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Abstract 

We compute the hyperfine splitting rrij/^ — m Vc on the lattice, using both the 
Wilson and 0(a)-improved (clover) actions for quenched quarks. The computations are 
performed on a 24 3 x 48 lattice at (3 = 6.2, using the same set of 18 gluon configurations 
for both fermion actions. We find that the splitting is 1.83 jj times larger with the 
clover action than with the Wilson action, demonstrating the sensitivity of the spin- 
splitting to the magnetic moment term which is present in the clover action. However, 
even with the clover action the result is less than half of the physical mass-splitting. We 
also compute the decay constants f Vc and fj/^, both of which are considerably larger 
when computed using the clover action than with the Wilson action. For example for 
the ratio fJ/^/fp 1 we find 0.32 ^ J; with the Wilson action and 0.48 ± 3 with the 
clover action (the physical value is 0.44(2)). 



The Hyperfine Splitting Lattice computations of the vector-pseudoscalar mass-splitt- 
ings m D * — rriD and mj/^ — m Vc , using the standard Wilson action for the quarks in the 
quenched approximation, give results which are much too small Q, |. At (3 = 6.2, for which 
the inverse lattice spacing (a" 1 ) is approximately 2.7 GeV, the discrepancy is about a factor 
of 2 for mo* — mo and a factor of about 4 for mj/^ — m^. In this letter we compute the 
hyperfine splitting mj/^ — m ric using both the Wilson fermion action, 

- q(x + fl)(^ + r)Ul(x)q(x)]j (1) 

and the nearest-neighbour 0(a)-improved (or "clover") fermion action M, 

Sp = Sf - irgo^a 4 ^ q{x)F flu (x)a flu q(x) (2) 

with the same set of 18 gluon configurations in each case. The computations are performed 
on a 24 3 x 48 lattice at (3 — 6.2 with r = 1. These 18 configurations have been used 
earlier in our study of light quark spectroscopy and meson decay constants, the results and 
computational details can be found in ref . ML [S] . For the quantities studied in M the results 
obtained with the two actions were broadly compatible. However, here we show that the 
hyperfine splitting in charmonium is almost a factor of 2 larger with the clover action than 
with the Wilson action (see eq.(|3|) below), demonstrating the sensitivity of this quantity to 
the magnetic moment term in eq.(^). We also present the results for the decay constants of 
the rjc and J/ip mesons. 

Qualitatively similar results for the hyperfine splitting were obtained by the Fermilab group 
Q, who have performed simulations using a fermion action which is similar to that in 
equation but with a factor of 1.4 multiplying the second term. This factor is their estimate 
of the effects of higher order perturbative corrections, and was obtained using a mean field 
theory calculation 0. These authors also find a larger value of the hyperfine splitting 
with their action than with the Wilson action. However this comparison is obtained from 
computations on lattices of different size. Below we will compare our results with those of 
ref. I. 

For this study we took k = 0.1350 for the Wilson action and k = 0.1290 for the clover 
action. These values were chosen so that the pseudoscalar meson masses are almost equal 
for both actions, and correspond approximately to the physical mass of the rj c . In table [TJ we 
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Wilson 
clover 


1.066+ I 
1.071+ I 


1.076+ i 

— 5 

1.088+ I 


0.0104+ 8 
0.0190 + }2 



Table 1: Masses (in lattice units) of the pseudoscalar and vector heavy-heavy mesons. 



present our results for the masses of the vector and pseudoscalar mesons (in lattice units), 
and for their difference. Setting the scale from the string tension gives a -1 = 2.73(5) GeV 
H Q and using this value we see that the masses of the mesons are within a few percent 
of their physical values. The masses were obtained by fitting the correlation functions in 
the time range t = 14 — 20, and the mass differences were obtained by fitting the ratio 
of the vector and pseudoscalar propagators over the same range. The fits were performed 
taking the correlations between the values at different time slices into account, and the 
values of x 2 /d.o.f were acceptable. The errors presented in this letter were obtained using 
the bootstrap procedure, described in detail in ref. ||, ||. 

The main result of this letter comes from the entries in the third column of table [l], from 
which it is clear that the hyperfine splittings are very different for the two actions. We 
stress that the results were obtained using the same gluon configurations and with the same 
analysis techniques. It is therefore likely that a number of the systematic errors would cancel 
in the ratio, for which we find: 



\ clover 



— m 



lie I 



\ Wilson 



1.831(1 



(3) 



In fig. |T] we plot the values of my — m 2 P (where V and P represent vector and pseudoscalar 
respectively) as a function of the square of the pseudoscalar mass. We include not only 
the values for charmonium obtained from table [l], but also those for mesons composed of 
light quarks for three different light quark masses 0. We see that for small masses the two 
actions give similar results, but as the mass is increased a gap gradually opens, with the 
clover action giving a larger value for the hyperfine splitting. 

It is interesting to note that the suggestion that this quantity might be susceptible to lattice 
artefacts has been made previously by the APE collaboration in a comparative study pf of 



1 The values of the inverse lattice spacing obtained by comparing the lattice values of the masses of the 
light hadrons to the physical ones lie within 15% of the result from the string tension. We take this as an 
indication of the size of the systematic uncertainty. 
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Figure 1: my — m 2 P versus mf>, in lattice units, for the Wilson and clover actions. 



Wilson and staggered fermion actions, albeit at stronger coupling and lighter quark masses 
than the present work. 

The result for the hyperfine splitting of charmonium obtained with the clover action is still 
much smaller than the experimental value. Taking a" 1 = 2.73 GeV, the values in table p] 
correspond to: 



28 + \ MeV 
52+ \ MeV 



Wilson Action 
Clover Action 



(4) 
(5) 



to be compared to the experimental value of 117(2) MeV. The errors quoted in (^) and (|5|) 
are statistical only, and the reader should bear in mind the uncertainty in the value of the 
lattice spacing. The corresponding values found by El-Khadra et al. are: 51(3) MeV at 
13 = 5.7, 62(4) MeV at (3 = 5.9 and 68(5) MeV at /3 = 6.1. Extrapolating linearly in a 2 to 
the continuum limit, these authors quote: 



mj/ip — rn Vc = 73 ± 10 MeV 



(6) 



where the error includes an estimate of the systematic uncertainty. The lattice spacing in 
this work was determined from the IP-IS mass splitting (a quantity which is considerably 
less sensitive to the form of the action 0). Thus it appears that there is a difference of 
about 20 MeV due to the different action used in ref.fHl. 
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Decay Constants The (dimensionless) decay constant of the J/^-meson is defined by: 

2 

(0|c(0) 7M c(0)|J/^) = e^ (7) 

where e M is the polarisation vector of the J/ip. The measured value of the decay constant 
is = 0.124(5). In our computations we take for the lattice vector current, the local 

operator Zy / c(0)7 /J c(0) when using the Wilson action, and the "improved" current 

Z% c(x)(l + y 7 - 5)7^(1 - yT D)c(x) (8) 

when using the clover action. Zy and Zy are the renormalisation constants (which ensure 
that the currents are correctly normalised), and can be evaluated in perturbation theory. 
We obtain the results 

"L 1 = 0.152+ \ Wilson Action (9) 

Zy Jj/il> 

-=^7— = 0.179+ 8 Clover Action (10) 

Zy fj/i> 

The Zy's have been calculated to one-loop order, Z|f=0.83 and 2^=0.90 if the lattice bare 
coupling constant g$ is used as the expansion parameter of perturbation theory, whereas 
Zy =0.71 and ^^=0.83 if an "effective" coupling g 2 e ^ = 1.75<7q is used (following suggestions 
in ref.0). For example using the effective coupling we find 1/ fj/f = 0.108 + | using the 

Wilson action, and 1/fj/j/, = 0.149 + I with the clover action. Hence the value of 1/fj/j/, is 
about 30-40% larger with the clover action than the Wilson action (although the uncertainty 
in the values of the renormalisation constants should be borne in mind). 

The decay constant of the r) c has not been measured, however we present the lattice results 
in order to compare the values obtained with the two actions. The decay constant is defined 
by 

1(015(0)7^750(0)1^))! = /^ (11) 

(with such a normalisation f m ~ 132 MeV). For the lattice axial current we take the same 
operators as in the vector case with 7^ — > 7^75. The corresponding results are: 

' f Vc a = 0.130+ % Wilson Action (12) 



z wJVo- ~ — _ 5 
-=cfn c a = 0.149 + § Clover Action (13) 

Z A 
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and Z% are equal to 0.87 (0.78) and 0.98 (0.97) when the bare (effective) coupling is 
used as the expansion parameter. Thus the values of f Vc are also about 30-40% larger when 
determined using the clover action than those obtained with the Wilson action. 

From eqs.(|9|)-(113f) we note that both the quantities l/Zy fjL and I/Za/tjc are larger when 
obtained using the clover action than the Wilson action. This is opposite to the results for 

the corresponding quantities for the light mesons tt and p [|J. Thus the difference obtained 

/-i 

for the decay constants with the two actions is amplified significantly in the ratios -pr and 

jp 

Using the chiral extrapolations for /„. and J" 1 , we find 





= 0.32 + 
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= 0.48 + 
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jp 
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where the physical value of this ratio is 0.44(2), and 

fvc 



■> Ve — o q + 5 



4.0 + 1| 



Wilson Action 
Clover Action 



Wilson Action 
Clover Action 



(14) 
(15) 



(16) 
(17) 



In these ratios the dependence on the renormalisation constants cancels. The differences 
between the results for the two actions in eqs. (|l4|) -(|T7|) indicate significant errors due to the 
fmiteness of the lattice spacing (at least, presumably, for the Wilson action) in the decay 
constants for the charmonium system. 



Conclusions In this letter we have shown that the value of the hyperfme splitting in 
charmonium in lattice simulations is very sensitive to the fermion action which is used, and 
in particular to the magnetic moment term in the improved action. Using the same gluon 
configurations, we have found that at (3 = 6.2, the ratio of the splittings for the clover and 
Wilson actions is about 1.8 (see eq.(|3|) ). Even using the clover action, the value we obtain 
for the hyperfme splitting is only about one half of the physical value. The clover action is 
a "tree-level improved action", i.e. there are no errors of 0(a), but the leading remaining 
errors due to the finiteness of the lattice spacing are of 0(a s a). Presumably at least some 
of the discrepancy between the value in eq.([|) and the physical one of 117(2) MeV is due to 
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these remaining 0(a s a) errors. El-Khadra et al. have tried to reduce these by performing 
a mean field calculation to estimate the effects of the higher-order perturbative terms on 
the magnetic moment term in the action || |7|. Their result for the hyperfine splitting, of 
about 73 MeV, although larger than that in eq.@, is still considerably smaller than the 
physical value. In view of the sensitivity of the splitting to the magnetic moment term in 
the action, it is likely that at least part of the discrepancy is still due to the finiteness of the 
lattice spacing. Unfortunately it is not possible at present to determine how much of the 
discrepency is due to the inadequacy of the mean field calculation, and how much to other 
systematic errors (such as quenching). We are forced to accept that the hyperfine splitting is 
currently not amenable to an accurate lattice determination. Nevertheless it is reassuring to 
find that one of the few quantities for which lattice computations with the Wilson fermion 
action give a result which disagrees significantly with experiment, is unusually sensitive to 
known systematic errors. 

For the decay constants of charmonium we also found significant differences between the 
values obtained using the two actions. For the ratio fj/lp/fp 1 (see eqs.([P|) and flT5] ) ) we 
find a result which is about 75% of the physical one with the Wilson action (consistent with 
the results in ref.0), and a result which is consistent with the physical one with the clover 
action. For the decay constant of the r\ c (see eqs. ([H]) and flTTD ) we find a larger result with 
the clover than with the Wilson action (note however that the error for the clover action is 
large, this error is dominated by the uncertainty in the extrapolated value of /„-). 

It will be very interesting to repeat this study for heavy-light mesons. In particular simu- 



lations with Wilson fermions indicate a substantial violation of the scaling law fp\/Mp ~ 
constant (up to logarithmic corrections) || [K| for the decay constants of heavy-light pseu- 
doscalar mesons P. This would imply a value of /# of about 200 MeV (one which is consistent 



with the simulations using the static approximation JTT], [T2|]), larger than many earlier expec- 



tations. It is now important to check whether these results will be stable under the reduction 
of the errors of 0(a) achieved by the use of the clover action. 



Acknowledgements This research is supported by the UK Science and Engineering Re- 
search Council under grant GR/G 32779, by the University of Edinburgh and by Meiko 
Limited. We are grateful to Edinburgh University Computing Service for use of its DAP 608 
for some of the analysis and, in particular, to Mike Brown for his tireless efforts in main- 
taining service on the Meiko i860 Computing Surface. CTS acknowledges the support of the 
UK Science and Engineering Research Council through the award of a Senior Fellowship. 



6 



References 

[1] M. Bochicchio et al, Nucl. Phys. B372 (1992) 403 

[2] A. Abada et al, Nucl. Phys. B376 (1992) 172 

[3] B. Sheikholeslami and R. Wohlert, Nucl. Phys. B259 (1985) 572 

[4] The UKQCD Collaboration: C.R. Allton et al., Phys. Lett. B284 (1992) 377 

[5] The UKQCD Collaboration: C.R. Allton et al., 0(a)-Improvement for Quenched Light 
Hadrons at (3 — 6.2: Mass Spectrum and Decay Constants, in preparation. 

[6] A.X. El-Khadra, "Charmonium with Improved Wilson Fermions II: the Spectrum", 
Fermilab Preprint FERMILAB-CONF-92/10-T (1992) (to be published in the pro- 
ceedings of Lattice 1991) 

[7] A.X. El-Khadra, G. Hockney, A.S. Kronfeld and P.B. Mackenzie, "A Determination of 
the Strong Coupling Constant from the Charmonium Spectrum", Fermilab Preprint 
FERMILAB-PUB-91/354-T (1991) 

[8] The APE Collaboaration: S. Cabasino et al, Phys. Lett. B258 (1991) 195. 

[9] CP Lepage and P.B. Mackenzie, Nucl.Phys. B(Proc.Suppl.)20 (1991) 173 

[10] C.R.Allton et al., Nucl. Phys. B(Proc.Suppl.)20 (1991) 504 

[11] C.R.Allton et al., Nucl. Phys. B349 (1991) 598 

[12] C.Alexandrou et al., Phys. Lett. B256 (1991) 60 



7 



